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ABSTRACT 


The use of composite laminates as structural members in various industries like 
aerospace, nuclear and automobile, has enormously increased due to their lightweight, 
high strength, high stiffhess, excellent fatigue resistance and better thermal expansion 
coefficients. All engineering materials inherently show variations in material properties. 
Composite materials experience larger dispersions in their material properties compared 
to conventional (isotropic) materials. This is because of more number of parameters 
involved in their manufacturing/fabrication process. Change in factors such as fiber 
orientation, laminate thickness, volume fraction, curing temperature and pressure, voids, 
impurities, curing time etc. induce variations in the lamina properties. This leads to 
variation in stiffhess coefficients of the laminates and thus results in uncertainty in the 
response behavior. The effect of variation can be accommodated in the formulation by 
modeling the material properties as random. The external loading is subject to variations 
due to uncertainties in the environment and is random in nature. Therefore, these 
variations in the stiffhess parameters and external loading have to be identified in the 
probabilistic sense for accurate analysis of composite laminates. The influence of 
randomness in material properties and loading on buckling, natural frequencies and 
deflection behavior of laminated composite shells has been analyzed in the present study 
with perturbation technique. The basic formulation of the problem is developed based on 
Classical Laminate Theory. 
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CHAPTER 1 


INTRODUCTION 

1.1 GENERAL INTRODUCTION 

Composite materials are extensively used by various engineering industries in 
various structural applications. Now a days, aircrafts are built with a very high 
percentage of components made from composite materials. Composite materials have 
many advantages over conventional materials, because of their superior performance in 
terms of their strength to stiffriess ratio. Composite laminates have the added advantage 
that their structural properties may be tailored according to the design requirements. This 
creates new possibilities and challenges for the analyst and the designer and makes 
accurate analysis a necessity for sensitive applications. 

The mechanical and physical properties of these materials are dependent on the 
production and fabrication conditions. All materials have dispersion in material 
properties due to lack of complete control over manufacturing/fabrication/processing 
techniques generally employed. Composite materials experience larger dispersions in 
their material properties compared to conventional (isotropic) materials due to larger 
number of parameters involved in their fabrication. In addition, the uncertainties 
involved in manufacturing and processing techniques employed for composites are more 
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as compared to techniques used for isotropic materials. These uncertainties can occur 
due to air entrapment, delamination and lack of resin, incomplete curing of resin, excess 
resin between layers etc. 

The conventional methods of design and anlysis of structures assumes the various 
parameters to be a constant. This approach ignores the variations in the system 
parameters and may not be acceptable for sensitive applications. The variations in the 
lamina properties result in variations in the stiffness coefficients of the laminates. 
Therefore, these variations in the stiffness parameters have to be identified in the 
probabilistic sense for accurate analysis of composite laminates. 

The randomness occurs mainly in the following two groups of parameters: 

• Geometric parameters like fibre orientation ‘ 0 ’, and lamina thickness ‘th’etc; 

• Material property randomness like elastic modulus, Poisson’s ratio, etc. of both 

fibre as well as the matrix. 

Considering the above aspects, the properties of structural materials may be modelled as 
random variables. 

In general, we may not have enough data about the conditions of the external 
loadings in engineering problems. In such cases, the exact estimation of these forces may 
not be feasible and it may only be possible to describe these as random. Therefore, the 
problem formulation with random external loading is more general, since it attempts to 
describe the system behaviour for whole class of possibilities'rather than for a single one. 

When geometric or material properties and external loading are random in nahue, 
the derived response parameters like deflections, natural frequencies etc. are also random, 
being functions of these basic system parameters and random loading. Depending on the 
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characteristics of the input random variables, the statistics of the response parameter can 


vary. 

In such conditions, statistical analysis not only provides an accurate method of 
analysis of the structural behaviour and determination of probability of failure for a given 
structure but also it can provide useful information on how to modify the structure to 
improve its performance. 

SHELLS 

The salient features of shells, as compared to other stmctural forms such as beams, 
frames and plates can be 

• Efficiency of load carrying behaviour; 

• High degree of reserved strength and structural integrity; 

• High strength to weight ratio; 

• Very small thickness ratio to other dimensions like span and radius of curvature; 

• Very high stiffiiess; 

• Containment of space. 

The usage of fibre reinforced and laminated composite materials have 
enormously increased in the domain of application and range of structural efficiency of 
shell forms. Ground, as well as space vehicles, having shell forms, have been designed 
and successfully built of high strength temperature resistant composite materials. The 
skin of aircraft structures and ship hulls are composed of shell forms, built of stiffened 
shells with composite material bodies. 
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SANDWICH CONSTRUCTION 


Instead of increasing the laminate thickness, sandwich construction is frequently 
used in aerospace applications. Sandwich with aluminium honeycomb cores or flex cores 
are used extensively in aerospace, land and water transportation and sports industries 
because of their high specific strength and stiffness. Some examples of strucfiires 
composed of sandwich panels are launch vehicles and satelliteshells, helicopter rotor 
blades, aircraft control surfaces, ship hulls, jet engine nacelles and skin, high-speed 
trains, etc. 

This type of construction consists of two thin facing layers separated by a core 
material. Several types of core shapes and core materials have been applied to the 
construction of sandwich structures. The core layer is made of low specific weight 
material like balsa, porous rubber, corrugated metal sheet, metallic and non-metallic 
honeycomb, etc., which may be much less stiff and strong than the face sheets. 

In addition to the possibility of achieving high flexural-stiffiiess providing a 
smoother aerodynamic surface in a high-speed range, sandwich-type constructions also 
exhibit many properties of exceptional importance for aerospace and civil constructions. 
Among these are 

• Excellent thermal and soimd insulation; 

• A longer time of exploitation as compared to stiffened-reinforced structures; 

• Possibility of being designed as to meet very close thermal distortion tolerances 
such as those required for communication satellite antennas and reflectors. 
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The most popular core is the honeycomb construction that consists of very thin foils in 
the form of hexagonal cells perpendicular to the facings. Face sheets usually consist of 
aluminium or fibre-reinforced composite laminates. Face sheets are stiff and strong 
because they carry most of the loads, while lightweight core separates face sheets so that 
a higher bending stiffiiess of the composite panel can be achieved. 

Regardless of the face sheet and core materials, the sandwich panel is considered 
a composite stmcture because of its inherent inhomogeneous and anisotropic nature. 
Even though the concept of sandwich construction is not very new, it has only been 
adopted for primary stmcture parts very recently. This is because there are a variety of 
problem areas to be investigated and overcome when the sandwich constmction is 
applied to design of dynamically loaded stmctures. To enhance the attractiveness of 
sandwich constmction, it is thus essential to have an understanding of the local strength 
and dynamic characteristics of individual sandwich panel/beam members. 

Advanced sandwich-type constmctions imply the presence of thick orthotropic 
core with bonded anisotropic face sheets that are treated as composite laminates. This 
arrangement presents an opportunity to tailor both the physical and mechanical properties 
of the faces by proper selection of laminate materials, their stacking sequence and fibre 
orientation. Suitable selection of fibre orientation and stacking sequence can result in 
substantial improvements of the buckling strength and the response behaviour to a variety 
of load conditions. The transverse shear elastic module of the core layer can also be 
optimised to enhance the overall response behaviour of the sandwich constmctions. As 
expected, analytical modelling of sandwich-type panels is much more intricate than that 
of the usual laminated composite stmcture. In contrast to the regular laminated 
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composite structure, for which the assumptions are postulated for the structure as a 
whole, in the case of sandwich-type constructions the asstnnptions involving the core 
layer are different from those associated with the face sheets. More over, the analysis of 
sandwich panels featuring laminated face sheets is much more complicated than that with 
single layered faces. The complexity is due to the presence of three types of asymmetries 
resulting from the lay-up sequences in the face sheets, namely: 

• Asymmetry with respect to the mid-surface of the face sheets, referred to 
as face asymmetry, inducing face bending-stretching coupling; 

• Asymmetry with respect to the mid-surface of the core, refer to as global 
asymmetry, which induces global bending-stretching coupling; 

• Presence of ply-angles between the principle axis of orthotropy of the face 
sheets materials and the geometric axes of the panel, inducing a structural 
coupling between stretching and shearing. 

1.2 LITERATURE SURVEY 

Any structural analysis problem has combinations of the following three major 
categories. 

1) Analysis involving different types of material; 

2) Analysis involving different types of structures shape wise; 

3) Analysis involving different types of loading. 

These lead to three types of problems with probabilistic modelling - with 

randomness in material properties, geometry and external loading. 
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Nigam and Narayanan [1], Crandall [2], Lin [3] have studied various class of 
problems related to external loading as random while the system properties and geometry 
are deterministic. 

Ibrahim [4] has presented a review of structural dynamics problems with 
parameter uncertainties. Manohar and Ibrahim [5] have reviewed number of topics on 
various formulations and solution techniques for structural dynamic problems with 
parameter uncertainties beyond the work reported in Ibrahim [4]. Vaicities [6] has 
obtained results for free vibrational analysis of beams with mass and flexural rigidity as 
random variables. Chen and Saroka [7] have investigated the response of a multi-degree 
dynamic system with statistical properties to deterministic excitations. They have used 
the perturbation technique to solve the governing equations of motion. Vinson and 
Sierakowski [8] have explained the behaviour of structures composed of composite 
materials. Leissa and Martin [9] have studied the free vibrations and buckling of the 
composite plates with variable fibre spacing. Zhang and Chen [10] have analysed 
complex stochastic structures subjected to arbitrary deterministic excitation. Salim et al. 
[11] have applied the perturbation technique for the static analysis of composite plates 
with various boundary conditions with elastic modulus, poison’s ratio as basic random 
variables. Also, Salim et al. [12] have investigated the effect of basic random parameters 
on the natural frequencies of rectangular composite plates using the perturbation 
technique. Singh et al. [13] have studied the initial buckling, natural frequency of 
cylindrical panel and composite plate with random material properties and obtained the 
second order statistics of response. Astill et al. [14] have examined the problem of impact 
loading of structures with random geometric and material properties. Elishakoff et al. 
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[15] have developed deterministic governing equations and boundary conditions for mean 
and covariance functions of displacement for statically determinate beam with 
stochastically varying stiffhess but with deterministic loading. Adah et al. [16] have 
discussed the problem of maximization of critical buckling load of an angle-ply 
laminated plate, when he properties of the plate are known to be scattered about some 
mean value. Yadav and Verma [17,18] have investigated free vibration for cylindrical 
shells with random material properties. Bucher and Brenner [19] have presented a 
method for analyses of stmctural systems with random properties subjected to dynamic 
loading. Low [20] has developed a fundamental but reliable and comprehensive approach 
of obtaining the exact roots for the frequency equation of beam systems. 

Jagadish [21] has studied the static response of graphite/epoxy and glass/epoxy 
composite plates with random material properties to random loading using the classical 
laminate theory. Onkar [22] has analysed the non-linear response behaviour of composite 
laminated plates with random material properties and random external loading. 

Gibson and Ashby [23] have explained the behaviour of cellular solids stmcture 
and their properties. Thomson [24] has presented a high-order theory for the analysis of 
multi-layer sandwich panels. Paik et al. [25] have analysed the strength characteristic of 
aluminium sandwich panels with aluminium honeycomb core theoretically and 
experimentally. Muc and Zuchara [26] have presented the analysis of a thin-walled 
sandwich laminate composite face subject to axial compression. Becker [27] has 
analysed, in closed-form, the thickness effect of regular honeycomb core material. Glass 
et al [28] have analysed Graphite/epoxy honeycomb core sandwich permeability under 
mechanical loads. Marimuthu [29] has studied the free vibration analysis of cylindrical 
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honeycomb sandwich laminate with random material properties using finite element 
method. Vinson [34] has explained the behaviour of shells composed of isotropic and 
composite materials. 

1.3 PRESENT WORK 

The behaviour of composite laminates with both kinds of randomness, namely, 
random material properties and random external loading, has been studied for a limited 
class of problems. The sandwich composite construction with random material property 
and loading has not been fully addressed. In aerospace structures, such situation occurs 
very frequently in many applications. Some examples are aerodynamic, jet and rocket 
exhaust noise, acoustic loading on aircraft, rocket and spacecraft; panels. The water 
pressure on ship and submarine panels during high-speed operation and at the event of 
underwater explosion also experience loadings that are random in nature. The acoustic 
pressure loading due to explosion and launching of rockets and firing of gims on vehicle 
and building panels are also random. 

The aim of present study is to investigate the stability and dynamics of honeycomb-cored 
composite shells of various shapes like cylinder, cone and sphere and finding the second 
order response statistics. The material properties and external loading are treated as 
random in nature. In the present work, the stiflfiiess properties of lamina such as Eii, E 22 , 

LJjj , are considered as the basic material variables, and are assumed to be random. The 
study is used to infer the performance of a typical satellite payload fairing. 
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CHAPTER 2 


FORMULATION 

2.1 INTRODUCTION 

This section presents a general approach for analysis of composite laminates involving 
random material properties and random loading. For the formulation of the problem the 
basic material properties like longitudinal modulus, transverse modulus, Poisson’s ratio 
and the external loading are taken as random in nature and all other parameters are 
treated as deterministic. 

2.2 CURVED COMPOSITE PANELS 

Figure2.1a and 2.1b show the geometry along with the stress resultants for a shell 
element. Let (x, y, z) denote the orthogonal curvilinear coordinates (or shell coordinates) 
such that the x and y curves are lines of curvature on the mid-surface z=0, and z curves 
are straight lines perpendicular to the mid-surface. Ri and Ra denote the values of the 
principal radii of curvatures of the mid-surface. The lines of principal curvature coincide 
with the coordinate lines. 

The shell under consideration is composed of a finite number of N orthotropic 
layers, zk and zk-i be the top and bottom z coordinates of the k**" lamina. 
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2.2.1 THE DISPLACEMENT FIELD MODEL: 


The displacement field relations for composite curved panels under consideration are [30] 

— 2 3 

m(x, y, z, t)= (1+z/ri) u +z. 1+ z q)\+z 6 \ 

- 2 3 

V (x, y, z, t) = (l+z/n) v +z. ^ 2 + z q) 2+z 9 2 

wCx, y, z, t)= w 

( 2 . 1 ) 

where ‘t’ is time, w, vand w are displacements along the x, y and z coordinates 
respectively. u,v and w are the displacements of a point on the middle surface and ^ 1 
and ^ 2 are the rotations at z=0 of normal to the mid-surface with respect to the y and x 
axes, respectively. Coefficients (p\, q) 2 , 6 \ and 62 have system dependent values. 

Following the development by Shankara and Iyengar [31] for composite plates, a 
procedure is outlined here for curved panels. The displacement field is expressed using 
fiinctions as coefficients 


( 2 . 2 ) 


(2.3) 


u = (l+z/ri) u +fi (z). <l)\ + h (z) 

V = (l+z/n) V + fi (z). <p 2 + f 2 (z) 

w = w 

The above equations can be written as 


dw 

dx 

dw 

dy 


u 


(1+z/ri) u +fi (z). ^ 1 + fi (z) 0 1 


V = (l+z/n) V + fi (z). ^ 2 + fi (z) ^ 2 


Vt; = W 


where 9 1 = — and 92 ^ ^ 


dx 


dy 
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The displacement vector for the model is, A = ^,>^ 1 ) 

T 

For Classical Laminate Theory, the displacement vector becomes A = («, v, w) (2.4) 

2.2.2 STRAIN DISPLACEMENT RELATIONS 

The strain displacement relations using equation (2.2) referred to Cartesian 
Coordinate system is 


£■, = £° + z{k° +z^k^) 


£, =€° +Z^k!, 


£. = S° +z(k^ +Z^K2) 


£•5 = £-5 + z V] 


f g + z{k^ + Z~K^ ) 


(2.5) 


0 du w 


where e. = — + — ; 
' dxr, 


Q dv w 
dy 


.0 - . 

^4 “ ^2 ’ 

dy 


dw 

dx 




0 _ ^ du 


^6 = 


’ dx' 


AT, 


dy 


2 . 


0 _ d^2 , . 

#v ^ "i™ « 

dx dy 


,,2 ^ M . 

' 3h^^dx dx^^' 


2 ^ -4 d^w 

' dy 


* 3h^ dx dxdy dy 


, I -4,, 9w. 


1 -4., 


( 2 . 6 ) 
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Equation (2.5) gives the strain-displacement relations corresponding to the displacement 
field with mid-plane strains and curvatures. 

2.2.3 STRESS-STRAIN RELATIONS 

The stress-strain relations for k*'' lamina in the material coordinate axis, whose fibers are 
oriented with respect to ‘x’ is given as [31] 




'Qu 

Qn 

0 

0 

0 ' 


r > 




Qn 

0 

0 

0 




> = 

sym 



0 

0 

•< 

















to 

1.^ 

1 




(2.7) 


where 7’ and 7’ denote the longitudinal and transverse directions of the lamina, 
cr, , (7, , s,, and s, are stresses and strains in the direction parallel and perpendicular 

to fiber direction. a^, a,., Si^ and are shear stresses and shear strains in the 
respective planes. 

The stiffness coefficients are defined in terms of material properties as [32] 

Q =_Al_ Q = — 

22 , =012 Q66=Gn Qm=G,, 

1^12 = ^21 

where and are longitudinal and transverse elastic moduli, is in-plane shear 
modulus, G ,3 , G ,3 are the out-if-plane shear moduli and ’s are the Poisson’s ratios. 


^12 ■“ 


(l-^12''2l) 


Qss=G 


23 
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The transformed stress-strain relations for the k* lamina with respect to the laminate 
coordinate system are [32] 


'<^l' 


'Qu Qn Qu 0 o' 

r 

£1 



-rK -rK 


«^2r 


Q 22 ^26 0 0 

£2 

(^6 

^ = 

1 

toi 

g: 5^ 

0 

0 

‘£5 > 



— K — K 


0’4 


Qm 045 

^4 

V. -4 J 


0"3. 



““ jtr 

where Qy are the material constants of the k* lamina coordinate system [31,39], 

2.2.4 EQUATIONS OF MOTION: 

The equations of motion for curved panels subjected to in-plane loads, Nx, Ny, Nxy 
and distributed transverse dynamic loading ‘q’ including effects of transverse shear and 
rotary inertia may be written using the principle of virtual work [33], 


dN^ dN^ — 7 dw 
dX dY ' ^ dX 


dN^ dN^ Ti ■■ Tf X 77 ^ 
dX dY ^ dY 


dQ, ,52, 4,dK, 4 


dx'^dY 


-r^C 


dX dY 


-)+r-7T( 


3h^ ax" dY 


^+2 ^)- 


axar 


2 .. 

+ ^;r(T;7T) + 


,a w. ,, . aV . -dii d<j), , div , d$^ 


R, /?2 ^'ax" 

16 ^a^w a^v. 


ar^ ^ dXdY dX dX dY dY ' 9h^ dX^ dY^ 


aM, aM« 4 4 aPj a?^ - 7 - aw 

ax ar ^ ^2 1 7 ^^2 er ^ 4^1/5 


aMx ai/. 


4 ,ap. ap, 




■ + ^ - 22 +7-7^2 - r-V(— + — ) = I? + />2 - 7^5 

ax ar a' " 3 A^ ax ar " 'ar 


( 2 . 10 ) 


where q (x, y, t) is the distributed transverse load and N, , etc are the stress resultants. 
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The Inertia’s /■ and /,• , i=l,2,3,4 and 5 are defined by [30] 


7 r 2/3 


/, 47. 


47, 


^ R, ^ 3h^ ^3h^R, 


h=i,+r 


Ih. 

R. 


r - 7+^3 47, 47, 

■' 2 =i2 +T — YTTr-y 


R^ 3^ 3h^R^ 


h = 


4/4 47, 


3h^ 3h^R. 


77 _ 47 47, 

1 3 = — ^ + - 


3h^ 3h^R, 


- , 87, 167, 

I A =7, -y — |- + y — p 
' 3/!^ ’ 9 / 2 " 


8 ^ ,.167, 


I'^=h-YTA + Y 


3h^ W 


- 47, 167, 

3/2^ 9 / 2 " 


(2.11) 


(7„7,/3,7„7„7,) = Eg-V^(l,Z,Z^Z^Z^Z®)dz 


/C=I 


where, p* is the mass per unit volume of the k"' lamina. 
2.2.5 THE LAMINATE CONSTITUTIVE EQUATIONS: 


It can be clearly seen from equation (2.9) that stresses in a composite laminate 
vary from layer to layer and hence a statically equivalent force and moment system is 
required for analysis. The stress resultants are expressed as (30) 


(V,,M,,7].)= i:g-'c7,^(l,Z,Z^)& 

K=\ ^ 

K=\ ^ 

K^l ^ 

N, = A,s;+Byk;+yEgk/ 
M,=B,s;+D,k;+rF,k; 
p,=E,s;+F,k;^H,k; 


i=l, 2,6 


( 2 . 12 ) 


ij=l,2,6 (2.13) 
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Where, R. and are the laminate stiffness expressed as 


r=l ^ 


(2.14) 


The equations for Classical Laminate Theory can be obtained by setting 7 = 0. 

Substituting all the above equations for curved panels in t 6 the constitutive equations and 
then in to" (2.10), we obtain equilibrium equations in terms of displacements. The 
equations of motion may now be written as 

[l]{A}={?} (2.15) 


where A. = (u,v,wy and q = (0,0,qy for CLT and L is a matrix of differential 
operators. The above equations are valid for laminated composite curved panels. 

2.3 CYLINDRICAL PANELS 

2.3.1 EQUILIBRIUM EQUATIONS: 


The general curved panel shown in fig 2.1a takes the form of a cylinder when Ri=oo and 
R 2 is constant. Using the geometrical definitions of cylindrical panels given in fig2.1a, 
equations ( 2 . 10 ) become 


dN,dN, 


dX dY 


= Lu 






dX dY 


dX dY 

N2 T" 

^ ^ ^ ^ ^ ^ arsr ~ ^ 


where the modified inertias I'., are given as /, = /, and 1 1 = /, 

(i,,h)=h!^-^pHxzydz 


(2.16) 
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2.3.2 STRAIN-DISPLACEMENT EQUATIONS: 


The Strain-Displacement relations are obtained by using equations (2.5) and (2.6) and 
cylindrical panel geometry definitions (fig 1.1a); 


^5 = 


du 

dx 

dw 

~Ec 


K°=0 


0 dv w 

S 2 =- + - 


^ ^2 

dv du 


^6 = — + 


dx ^ 


fc°=0 


*4 




<=0 




sy dx^ 


' y % 




3y dy^ 


' y^dx^ 


^^6 = 


•4 .2d^w 


Sh dxdy 


(2.17) 


2.3.3 CONSTITUTIVE EQUATIONS: 

The constitutive equations are obtained by using equations (2.13) 

= m,=o (2, = o y=i,2,6 


(2.18) 


2.4 SPHERICAL PANELS: 

2.4.1 EQUILIBRIUM EQUATIONS: 


The general curved panel shown in fig 2.1a takes the form of a sphere when Ri= R2=R. 
Substituting the geometry definitions for spherical panels in to equation (2.10), we get 


dN, dN, 

— L4. 6 

dX dY 


= I^u 


dN. dN, 

— 14- — 2. 

dX dY 


= /',v 
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^1 ^2 -^1 ^2 \T XT XT I \ T ■■ 


dX dY R R ^'dX^ ' dY^' ^'^dXdY 


where the modified inertias /,', I. are given as /, = /, and I[ = /, 




(2.19) 


2.4.2 STRAIN-DISPLACEMENT EQUATIONS: 

The Strain - Displacement relations are obtained by using equations (2.5) and (2.6) and 
spherical panel geometry definitions (Figure!. la). 


du w 


dv w 


dw 


s = + — 

‘ dx r, 

s. = — ^ + — 
dy 

1 

II 


0 dw 

Sr — 

dx 

0 dv du 
= — + — 
dx dy 



K° =0 

K°=0 

^6=0 




= — 
‘ 3h^ 

2a=w 

dxdy 

, _-4 dw 

■>F<S> 


(2.20) 


2.4.3 CONSTITUTIVE EQUATIONS: 

The constitutive equations are obtained by using equations (2.13) 

N.=A,jS; M, =0 Q,=0 ij=l,2,6 


( 2 . 21 ) 


2.5 CONICAL PANELS: 

2.5.1 EQUILIBRIUM EQUATIONS: 

The general curved panel shown in fig 2.1a takes the form of a cone by setting Rl= oo 
and R2 varying linearly along the x-axis gives a conical panel. 
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Substituting the geometry definitions for conical shells in to equation (2.10), we get 


dX 8Y ' 




dX BY 


iv 


dX dY R, 


,8^w 


,8^w 


8^w 


+ C-^) + ^xr = I.^-q 


8X 


8Y 


8X8Y 


where the modified inertias I'., I. are given as /, = /, and I[ = /, 


(A.A)=zgrv'('-2)* 


K=\ 


( 2 . 22 ) 


2.5.2 STRAIN-DISPLACEMENT EQUATIONS: 

The Strain - Displacement relations are obtained by using equations (2.5) and (2.6) and 
conical panel geometry definitions. 


n _5u 

S'* ~~ 

8x 

0 Qw 

6. = 

8x 

0 5v w 

g _ j 

dy 

a dv du 

e: = — + — 

8x dy 

dy 


^ o 

11 

O 

kI=Q 

kI=0 


„2 _ -4 

Zh^^8x-^ 

-4,ax 

3h^^8y^^ 

2 -4 

^6 = 7 

' 3h^ 

(23 V 

8x8y 


, -4 5w 


(2.23) 


2.5.3 CONSTITUTIVE EQUATIONS: 

The constitutive equations are obtained by using equations (2.13) 


N,=AyS° M, =0 e,=0 i,j=l,2,6 (2.24) 
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2.6 GENERAL NATURE OF THE PROBLEM 


The mathematical statement of a structural analysis problem usually involves three types 
of terms as follows. 

1 . Terms describing known system parameters. 

2 . Terms describing the known forcing function. 

3. Terms describing the unknown response vector. 

In general, the governing equations of any structural analysis problem can be 
expressed as: [i]{a}={ 9 } (2.25) 

Where, [z] -System matrix (m X n), {a} -R esponse vector (n X 1) and {^}-Loading 
vector (m X 1). 

In the above equation the characteristics of [z] and { 9 } are known while those 
of {a} have to be determined. [z]can be expressed as known functions of a set of 
primary variables ( 6 , ). The forcing vector is random in nature and can be 
expressed as known functions of set of loads . Elements of {a} are unknown and 

also are the functions of ( 6 / ) and J . The problem is therefore, is to find the statistics 
of {a} when the statistics of the primary random variables (p , ) and the statistics of the 
applied loads (g^ are known. 

2.7 RANDOM SYSTEM PROPERTIES AND LOADING 

The randonmess in material properties and loading are taken into account using 
‘perturbation technique’. It is assumed that all the primary variables and loading 


21 



cx)mponents are independent of each other. Also it is assumed that the dispersion 
(randomness) of each random variable about its mean is small in comparison to the mean 
value. This is common to most engineering problems. 

Any random variable can be expressed as a zero mean random part superimposed 
over a mean. This can be expressed as 

RP = RP‘'^RP^ (2.26) 

The random process equation [i]{a} = can be represented in above form. The basic 
random process as well as the loading can also be expressed in this manner. 

Thus we have 

[l]~[l‘]+[u\ 

{a}={a'}+{a' 

bi=b^ +b' (2.27) 


where b , = basic random variable. 

As mentioned earlier, for many practical engineering situations the random component 
{RP') is small compared to mean {RP^). Therefore the equation (2.27) can be written as 
follows 

[x,]=[V]+4vl 

h;}=K/K^h/} 

bi=bi+d>[ (2.28) 
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Here ‘ ‘ is a perturbation, small in magnitude. No generality is lost in using the 
same value of ‘ ’ in all the cases. The remaining factor in the random term can absorb 


varying levels of this term. 

Using Taylor’s series expansion, some of the terms of equation (2.28) can be expanded as 


54 

eL,; = IL^sb; 

' I db,‘‘ 

^ dA/ ^ 5A / 
sA / = Z + 1 — 

' raw" Ta^/ 


(2.29) 


As fib/ are small in magnitude, the second and higher order terms are neglected here. 
This assumption is justified, if the randomness is small compared to mean value and sets 
a limit to the range of applicability of the approach. 


Substituting equation (2.28) in equation (2.25) 
[Z" +£Z:^]{a" +fiA^)=fe" ^sq^} 


(2.30) 


Expanding equation (2.30) and comparing terms of each order of ‘ fi ’, we get 

[l"A" + fit," A" + fiL^A" + s^Uk]= +sq'] 

Collecting terms of different orders of , we get for zero* order and first order terms, 

[LY{^Y=kY 

fi’ > [lY {aY+[lY{aY = {qY (2.31) 


Substituting equation (2.29) in equation (2.31) 


I I,-/ 


ba; 

/ dbl 


dA 


j ^ '' 


i dqj^ 


<lk 


+ S 

i 


az,/ 

I— 

T db, 


-WU 


=:q. 


(2.32) 
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Since dispersion about mean is small, the derivatives of the random variables can be 
approximated by the derivatives of the mean values. Then equation (2.32) becomes 




5A f ^ dA f 
I dbf * dq,' * 


+ E 

j 


dL. 




=qi 


(2.33) 


dA dA 

Here the unknowns are — -r- and — 


db, dq, 


dA: 


To find — ^ , multiplying equation (2.33) by bj and taking expectation, we get 
db, 




. 5a; 


^ dL.: 

d ^ ij 


J I " ab/ ' “ J * ^ dq," “ JI ' db, 


is the cross correlation between x and y. 


d bib^ q.bu 


(2.34) 


Assuming that the primary variables and the loading components are independent random 
variables, the above equation becomes 


J db^ J db^ 

This is satisfied only when 


ai,/ 


(2.35) 


d^^i 


. dL. 

d ^ ij 


db 


= 0 


(2.36) 


dA:" 


Solving the above equation, we will get — ^ . 

db" 

dA j" r 

To find — ^ , multiplying equation (2.33) by and taking expectation, we get 

d^^i' 

EZA/ — 

y T dq," 


(2.37) 


Again with assumption that 's and q^ 's are independent of each other, the above 


dA:" 


" " -K. . =K, 

^ dq: 


equation becomes 


(2.38) 
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5A/ 

Solving the above equation, we get — . 
Therefore, the total response becomes 


Aj=A/ +f| 


SA; dA; 


db; 




The variance of the response can be evaluated as follows 

Var(A,)=4^-V)' 


Var(A^.)=-£^ 


dA; ,5A, 

ZV— ^ + 2:9* — T 

/ db; * dq" 


* J 


( 2 . 39 ) 


( 2 . 40 ) 


Standard deviation is the square root of the variance. The above approach is general in 
nature, and can be applicable to different problems in analysis of structures. 
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CHAPTER 3 


STABILITY AND DYNAMIC ANALYSIS 

3.1 INTRODUCTION 

This chapter presents the general approach for stability, free and forced vibration 
analyses of laminated sandwich composite shells with random material properties and 
external loading. 

3.2 BUCKLING OF SHELLS 

The increasing use of high strength and stiffriess materials in recent times for 
structural component design and the need to produce light weight and optimized 
structures have lead to widespread adoption of thin walled components in aerospace 
structures. The aerospace structures, generally made up of lightweight materials, 
experience complex loading during service life. Instability may occur in these structures, 
leads to buckling failure. To avoid this, evaluation of critical buckling load is an 
essential requirement for design of optimized structural components. A unified approach 
is discussed for buckling analyses of cylinder, sphere and conical shells that involve 
anisotropic, layered composite construction. 

This section outlines a probabilistic methodology for application of classical 
approach in conjunction with first order perturbation technique (FOPT) for evaluation of 
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second order statistics of critical buckling load (mean and SD) as a function of the basic 
random variables - Longitudinal modulus, transverse modulus and Poisson’s ratio. 

Classical thin plate theory has been adopted for the analysis. According to this 
theory, a linear element of the shell, extending through the wall thickness and normal to 
the middle surface in the unstrained state, upon the application of load satisfies the 
following conditions; 

• Undergoes at most a translation and a rotation with respect to the original 
coordinate system. 

• Remains normal to the deformed middle surface. 

• It is assumed that shell resists lateral and in-plane loads by bending, transverse 
shear by membrane action and not through block compression or tension of 
the shell wall in the thickness direction. 

• The linear element does not elongate nor contracts. 

• It remains straight upon the shell being loaded. 

• Saint Venant principle applies. 

The main assumptions in the classical laminate theory may be put as 

• Thickness of the panel is small compared to its other lateral dimensions. 

• Plate deflections are small compared to its thickness. 

• Normal stresses in the transverse direction can be neglected. 

For buckling of the shell under external loaTding, it is assumed that pre 
buckled deformation are not taken into account and ends of the cylindrical shell are 
supported by rings rigid in their planes but with no resistance to rotation or bending out 
of their plane. 
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3.2.1 STOCHASTIC CLASSICAL APPROACH 


For all edges are simply supported, exact Navier type solution is possible for curved 
panels. When exact solution is combined wth probabilistic method, it is possible to 
analyze the random characteristics equation arising from buckling of composite 
laminates. In this section, a detailed study is presented to obtain the second order 
statistics of the buckling loads for sandwich cylindrical, spherical, and conical panels 
using stochastic classical approach (SCA). 

Boundary conditions are; 

Referring to the figurel.la 

x=0,a v(0,y) = v(a,y) = 0 

w(0,y) = w(a,y) = 0 

Ni (0,y) = N 2 (a,y) = 0 

y=0,b u (x,0) = u(x,b) = 0 

w (x, 0) = w(x,b) = 0 

Ni(x,0) = N2 (x,0) = 0 (3.1) 

The displacements satisfying all the boundary conditions can be expressed as 


= sin 

IB=1 n=l ^ O 

= cos- 

a b 


m=l n=l 

K^,j)=SZ^™sin— sin 
a b 


rriTDc . nuy 

a b 

mTDc njw 

cos— ^ 

a b 

. mjDc . ww 




(3.2) 


Substituting the equation (3.2) into equation (2.25) results in a homogeneous system of 


equations 


=0 


(3.3) 


}=’ 
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where ‘k’ is a constant column vector, i, j=l,2 and 3 for CLT. For nontrivial solution of 
‘ k’ from equation (3.3), the determinant of coefficients should be zero: 

|M,| = 0 (3.4) 

Equation (3.4) can be written in expanded form as 


11 

«I2 

(2,3 

21 

^22 

«23 


*^31 ^32 *^33 ^cr 


(3.5) 


Where the expression for critical buckling load is given as 

^cr = N I a~) + N TT' I b''-) + N I ab) (3.6) 

In the above equation, ay’s are functions of system stiffness, constants “m” and “n”. 
Expanding the above determinant gives the expression for Ncr in terms of ay’s. 

Ncr=F(ay) (3.7) 

The stiffness elements ay are random in nature, being dependent on the system material 
properties, consequently, the buckling loads are also random. 

3.2.2 SECOND ORDER STATISTICS OF CRITICAL BUCKLING LOADS - 
PERTURBATION APPROACH 

A random variable can be split up as the sum of its mean and the zero mean random part 

(3.8) 

Similarly, ay can also be expressed as a.j = + Uy (3.9) 

On substitution of (3.8) and (3.9) in (3.5), yields 
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(3.10) 


d 

r 

d , r 

d 


«I 1 

-t-a,, 

a , 2 +a ,2 

«13 

+ <^13 

d 

r 

d , r 

d 


^21 

+ ^21 

^22 ^22 

^23 

+ (3,3 

d 

^31 

+ 031 '' 

d , d 

^32 +"^32 

d , r 

^33 +^33 • 



Expanding (3.10), collecting the same order of magnitude terms, and keeping only up to 
first-order terms, one obtains the following relation in the symbolic form: 

Zero*’' order: 

Nc/=F(ai/) (3.11) 

First order: 

Nc/=F(ai/,ai/,Ncr‘’) (3.12) 

Equation (3. 1 1) is a deterministic equation relating the mean quantities. The mean value 
of buckling loads can be obtained by solving the following expression 


^ 1 , 

^12 

^13 

^ 2 : 

^22 

*^23 

^31 

^32 

^33 



a,. 

^12 


i 

« 2 I 

^22 



(3.13) 


Equation (3.11) gives an expansion 


f.T r ffd d r, d d r, d d r, d d d d r , r\ d d r 

^cr “((^11 ^33 ^22 "^^22 ^33 ^11 "^^11 ^22 ^33 ”^^12 ^23 ^13 ^23 ^12 ^\3 ^23 

d d 


d d r 

“ T 

'11 


d d 
^22 
d d 


ry d d r r d^ ^ ry d d r r d^ ^ d d r r c /2 

’^^11 ^23 ^23 ^23 +2 aj 2 ^23 ^13 ^22 ^13 ^^22 ^13 ^13 ”’ ^33 ^12 ) 

• Nj (a,/ a,, ^ “ -2a, /a,, ^))/G 


Where G = -< 2 , 2 “^ ) 


(3.14) 


Equation (3.12) is the random equation, which is a function of aij**, ai[ and Ncr*’- For 
present case, ai/ and Nc/ are random because the material properties are random, as 


30 



discussed earlier. Let b/, hz denote the random material properties. The bi"^ can 

also be expressed as 

b,=b,'^ +b[ (3.15) 

Expanding the above expression in Taylor’s series, when bi*^ is small compared with its 

mean value, one can expand Lij"^ and Nc/ about bi‘*=l,2. . .m. keeping only the first-order 

terms, one obtains 


/=! 


/=! 


(3.16) 


Where ‘j’ denotes the partial differentiation with respect to bi ; and the derivatives are 
evaluated at bi*^. 


bi=Eu, b 2 =E 22 and b3=L>,2 

Substituting (3.16), in to (3.14) and simplifying, we get 
NJ = PE,; +QE,,' 

From basic definition of variance. 


Var{N„) = E 






>/ 



(3.17) 


= ttNj„Nj,,cov(b'X) 


This results 


vaiiNj) = p^ 




SD(Ncr>Sqrt(var(NcrO) 

where 


(3.18) 


(3.19) 

(3.20) 


31 



<233 _p d ■'■‘^22 <^33 ^ d ■'■^11 ^22 d + 2 ( 3,2 <^23 “ 

Cr 5£,, 9£,] 9£„ 9£^ 

^ d ^ d ^ ^ ^ ^ 

, ^ d d oaj^ « d oa^2 r. d d oa^2 ^ d d ocLi^ 

+ 2 aj 2 aj 3 ^ + 2^23 ^ 2a^2 ^33 ^ 4 2^23 ^11 ^ 

9£^ii o^ii 9£„ 95^1 j 

^ 00,3“' rf2 aa,,'' ^2 5^2/ ^ rf2 aajj'' ^ 

„ rf “23 d “l3 d “12 dl 

E„ a£:„ a^„ a£„ 

d da,,‘‘ , aa,/ d 5a, j'' 

1 ^ ■^^22 ^ ^^12 4 si 

dE., 


rs d 4 

-2a, 3 a22 — ^ 

a£. 


;,;r 

a£. 


'll 


9i?, 


'11 


“'ll 


-'ll 


^ 1 P ^ 4 4 9^2 

0 = -[ki ^33 — 

G dE. 

^ d d ^^23 ,0 ^ ^^12 r% d d ^ d d ^^7 

+ 2 ( 2^2 ^13 ^ +2^23 ( 2 j 3 . _ ^ “ 2 aj 2 <233 •—- 7 *'- 2(323 a,, — — 

9^22 5^22 5^22 5^2 

9a. 


<f ^^22 ^ ^ 9a. j j ^ dcLy^ r\ d d 9a. 

7T“^^22 ^33 T7rT"^^ii ^22 T7rT“^^^i2 ^23 "T^ 

^22 5 JS 22 ^^22 

d 

d 


22 


23 


^ ^^13 ^^11 ^^22 

.5 -_ w "35 ^ ^13 d ^12 

9£^22 


”2aj3 a22 ——^-^23 “ ^ 

9^22 dE22 


42 d^ 
12 i/i 


- 

22 


\j d ( d 9a 22 ( 

k, 

dE,, 


d 5 a,^ , ^ 


dEj '' a£. " 


-'22 


'22 


D ^ r f ^ ^ 9 a 22 </ rf 9 a j I 4 4 9 a 33 , </ d 9 a. 3 

•^~ 7^!-{^11 <^33 "I 7 "^^22 *^33 "I 7 "^^11 ^22 "I 7 ^^^12 ^23 T i 

G do,, do,, do,, do,, 

« d d 9 a 23 ^ d d 9 aj 2 d 4 ^^\2 r\ d 4 ^^21 

+ 2a^2 ^13 T 7'^2a23 aj3 ~ j — 2a^2 ^33 T 7 ““^^23 ^n T 7 

9^,2 5^2 5^12 9^2 


9a,, 


42 9a ,, 42 9a 22 4 ^ 9a 33 

^ 7 ”^13 T 7 ”'^12 "I 7 / 

9^2 9u,2 9Uj 


9 ^ ^^13 

- za,3 a22 ^ ^23 

9^2 9^2 

■XT d f d 9a,- 4 9a,, 

■-^cr {^11 *7 7 ”^^22 T j — 2a^2 p 

9t;,2 9(;,2 9^. 


'12 


'12 


do. 


'kii 


12 


-"12 


(3.21) 


3.3 FREEVIBRATIONS: 

The free vibration equation is obtained by setting the in-plane loads and excitation 
terms to zero in the governing equation (2.10). 
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For all edges are simply supported, exact Navier type solution is possible. When 
exact solution is combined with probabilistic method, it is possible to analyze the random 

generalized eigen value problem associated with free vibration of composite laminated 
panels with random material properties. 

In this section, a detailed study is presented to obtain the second order 
statistics of the natural frequencies for sandwich honeycombed cylindrical, spherical, and 
conical panels using stochastic classical approach. 


The displacements satisfying all the boundary conditions can be expressed as 

cos^sm^expifca) 


(3.22) 


a b 

v(x,y,t) = exp(iM) 


fn—\ /r=l 
00 00 


«=sl 
00 CO 


.v(x,y,r) = 2;E»'..sin^sin^fexp(,ra) 

. «=i /i=i Cl b 

Substrtuting the equation (3.23) into equation (2.25) results in a homogeneous 


; system of 


equations 


(3.24) 


^Mgkj=Q 

Where k is a constant column vector, i. j=l,2 and 3 for CLT. For nontrivial solution of 
‘ k’ ,In equation (3.24), the determinant of coefficients should be zero: 

KI=o 


(3.25) 


Equation (3.25) can be written in expanded form as 



^12 

a , 3 

to 

^22 

^23 

^3, 

^32 

^33 


The above equation can be expanded as 


(3.26) 
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+I^X-I, =0 


(3.27) 


where X = paP’ 


~ ^22 ^33 

I ^ — Cl^iO.22 + <^22^33 ^11^33 ~ ^12^21 ~ ^23^23 ~ ^13^3! 

■^3 ~ ^11^22^33 ” *^11^23^32 ^12*^21^33 ^12^23^31 ^13^21*^32 ” ^31^22^13 (3.28) 

To obtain the statistics of the natural frequency, the variable terms can be 
written as sum of the mean terms and the corresponding zero mean random terms. 

A = A" + 


Likewise, all other random quantities are defined as follows 



d 

. r 


d 

r 



, r 


” ^11 

+ ^n 

^22 

■“ ^22 

+ ^22 

«33 

= «33 

+ <333 


d 

, r 



. r 


d 

. r 

a, 2 

^12 

+ a^2 

a, 3 

= ^13 

+ ^I3 

^23 

= 

23 

+ ^23 


Substituting these in the above equation (3.27) 


+ 3 X‘‘X'^ + 3 X^h^)-{x‘‘" +X'^ + 2 X‘‘xla,/ +a„' + 02 / + 022 '' +^ 33 “' +0^/] 

. I 'id , or If d d . d r, r d , r r. d d , d r, r d , r r 

4 -\A +A j{ajj ^22 +^11 ^22 "^^11 ^22 “^^ll ^22 “^^11 ^33 *^^11 ^33 ^33 ”^^11 ^33 

. d d , dr, r d , r r d d dr r d r r 

^22 ^33 ^22 ^33 ^22 ^33 ^22 ^33 ^12 ^21 ^12 ^21 ^12 ^21 ^ ^12 ^21 


d d dr r d r r d d dr r d r r^ 

~^23 ^32 “^^23 ^32 ~ ^23 ^32 ~ ^23 ^32 ~ ^13 ®31 ~ ^\3 <^31 “^13 ^31 ~ ^13 ^ 3 \ J 


34 



■{(^11 ^22 ^33 “^^11 ^22 ^33 "^ ^11 ^22 ^33 “^^11 ^22 ^33 "^^11 ^22 ^33 "^^11 ^22 ^33 


“^^11 ^22 ^33 "^^11 ^22 ^33 ) (^11 ^23 ^32 ^23 ^32 “^^11 ^23 ^32 *^^11 ^23 ^32 


, d d r , d r r ^ r d r ^ r r r\ / ddd d r d 

“^^11 ^23 ^32 "^^11 ^23 ^32 ^23 ^32 “^^11 ^23 ^32 ) “ (^12 ^21 ^33 "^^12 ^21 ^33 


"^^12 ^21 ^33 "^^12 ^21 ^33 "^^12 ^21 ^33 “^^12 ^21 ^33 “^^12 ^21 ^33 “^^12 ^21 ^33 ) "^ 


/ddd, d r d . r d d r r d d d r d r r w 

(^12 ^23 ^31 “*'^12 ^23 ^31 "^^12 ^23 ^31 "*"^12 ^23 ^31 "^^12 ^23 ^31 **" ^12 ^23 ^31 )} 


(3.29) 


Taking expectation of above equation, we get 


A d^ 'id^f^^.^d d ^ n J / d d d d d d d d d d d d \ 

^ {^11 ^22 ”^^33 } ^ {^11 ^22 ^11 ^33 "^^22 ^33 ^12 ^21 "”^23 ^32 ^13 ^31 } 

c^d d d ddd ddd ddd ddd d d d ^ r\ 

-{«n «22 «33 -«n ^23 «32 -^12 ^21 ^33 + ^.2 «23 ^31 + « I 3 «21 ^32 -^31 «22 ^13 }=0 

(3.30) 

As, the expectation of all random terms is zero. The above equation is the characteristic 
equation for the mean part of the random natural frequencies. Subtracting the equation 

(3.30) from (3.29), and on simplifying, we get 

^ =-^{^ (^n +^22 +^33 )-^ (^1. ^22 +«'n «22 +^11 «33 + «33 +^22 ^33 

+ CI 22 CI 22 ~ fl(j2 £^21 ~ ^12 ^21 ^23 ^32 ~ ^23 ^32 ~ ^13 ^31 ~ ^13 ^31 ) } 


(3.31) 


where 


G = 3A‘'" -2A‘'(ai/ +a22‘' +^ 33 "') 


/ d d d d . d d d d d d d d \ 

+ (^n ^22 ■^'^11 ^33 "^^22 ^33 ~^12 ^21 “^23 ^32 ~^13 ^31 ) 
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F ddr d d r d d r d d r d d r d d r 

"■^11 ^33 ^22 ”^^22 ^33 ^11 ’^^11 ^22 ^33 ““^11 ^32 ^23 ^23 ^32 ^11 ““^11 ^23 ^32 


d d r d d r d d r , d d r , d d r , d d r 

’^12 ^33 ^21 “"^21 ^33 ^12 "“^12 ^21 ^33 “^^12 ^31 ^23 *^ ^23 ^31 ^12 “^^12 ^23 ^31 


ddr d d r d d r d d r d d r d d r 

■< 3]3 ^32 <^21 ^21 ^32 ^13 ”^^13 ^21 ^32 ~ ^31 ^13 ^22 ~ ^22 ^13 ^31 ~ ^31 ^22 ^13 


(3.32) 


aj can be expressed as a function of input random variables, through Taylor’s series 


expansion as 


da. 3 da. 

ay = = I— V*/' 

" I dbf M dbl’^ ‘ 


(3.33) 


Taking b^ = £■„ ,£>3 = £'22 ^3 = t^n 

Using above equation, and simplifying, we have 

A'=F[£,.']+ek']+«k1 

where 


(3.34) 


G 

dE,f dE,, dE,,^ 


^ 33 ^ I _ d.da,f da^; 

V d’^ d) 

SEjj 


_ ^ I I l‘‘ Ua << ^22 I ^33 

G dEji dEjj dE^^ a^j, a£:„ 

1 rf d^23 di ^^13 \i r rf da.y 


, d d OG, 

““ iZ' 


'11 

^ ^ ^ ^ ^^33 , O,. ^ ^^13 , o« ^^23 , ^ ^ 1 : 

j j T £Zjj ^22 // -^^19 ^23 ^ ^^12 ^13 A ^*^23 ^13 

a^;,/ dE,^ dE,f '' a.E:,/ ' dE 

d 


„ d d . - d d da^. 

“2^2^ ^33 j 2^23 Gjj 

a£',/ a£;„ 

d 


n /T tf2 ^<^33 

^ “12 rfjJ 

dE,, dE„ 


'11 

^^23 ^ 

I OE, ^“13 “22 ^ ^23 ^ 

oE,, oE,, oE,, 
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n = _Lr , ^^22 , ^^33 

O' 9^22 ^-^22 ^^22 


j + rf) 1(^11 ;r^ + Tr-7) + “22 (T^rT + T^rr) 
9£^5, 9£„ 9£„ 9£^„ 9£^„ 


is ^ O 

.oa,, da,. 


J vwjj ^^22 \ o/ « 

+ ^33 (T^ + T^)-2(a„ —^ + ^23 

9^22 ^^22 <5£22 


5^22 


d d ^^\\ d d ^^33 ^ d d ^^13 ^ d d ^^23 /\ d d ^^]2 

+ '^’2 ^33 -p rf ■‘■^n ^22 TtrT + 2<3,2 O^ ^ + 2 ( 2,2 ‘213 ^,3 ^ 

0 /S 22 dE^2 dE^ dE^2 


d d 5a,, 


^ a a '-'“12 ^ d d ^^23 ^ d 

- 2 a,, a33 ,/ ~2a,3 a,, — — j--2a,3 a 

0^22 5£22 


</ 5a, 3 fj2 dct,t 

13 “■>■> ^-^23 ^ 

" dE„‘‘ dEJ 


" 5 ^ 22 " 


" dES^ 


^ r -jd^ , 503, Sajj ,rf (/_ d da^y , 5 a 


i? = -[r ( 
G 


— ^ + -^) - r {(a.,“ 

au,^'' au,3‘' au,3^' " 


d 1^ ^ 

22 I ^^33 


^.aa,,'' , aa33‘' 


+ -^) + 033 " 

5^2 5U,3" aO,3‘'' 


at>,2 at2,3 


+‘'./(?4+m-)-2«!,.'^ 


ai;,/ 


5^12 , 3«23 , , 9a,3 . 5a33 

t-t+«23 ;^-t+^.3 t-t)p{«u «33 7-^ 

5^2 5y,3 ai;,3 So, 3 


d 5a,, 

^22 «33 -T^ 

au„ 


I ^ ° ^ ^ f ^ ^^33 .0 ^ </ ^^13 ^ d d ^ d d 

“^^22 ^33 ^ // ^22 ^ +2a,5 - 3- + 2ao <2,, :r + 2a„ a,, 

ov„ 


d d^_ 
^12 ^23 ^ d 


23 "*I3 2s d 

ay, 3 


o y, ^ ^ ^^12 o >v ^ ^ ^^23 fy d d 5a, 3 ^2 aa,, 

- 2 a , 2 0^33 T-T-2a23 «„ T^-2a,3 ^33 -2V-a23 

au,2 ay,3 ay, 3 ay, 3" 


7 ±± /7 

M3 ^ d ^ d - 

ay, 3 ay,3 


Taking the second order moment of equation of (3.34), we have 

e[a^ ] = P^E E,{\ + Q^E £ 33 ^' +R^E E,^^ + PQeIe^.^E^^' ]+ 

It is assumed that all the four zero mean random variables are independent of each other, 
so that cross-correlation terms vanish. 

Thus, the above equation takes the form; 

j+«“.Ek'i (3.35) 

By definition of the variance, we have 
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,2 


(3.36) 


FarCr') = P'of:,/' + Q'^aE^f + P VC/,/ 

SD(A^') = +<2Vp2/' +PVt/,/' (3.37) 

Thus, using the above equation SD{a)'^) is obtained. 

3.4 FORCED VIBRATIONS 

3.4.1 INTRODUCTION 

In general, vibration is undesirable for all aerospace structures. When a structure is 
subjected to external excitation, it is forced to vibrate at the same frequency as that of the 
excitation. These vibrations may be undesirable for the structure whose functioning may 
be disturbed or it may lead to its fatigue failure, if large vibration amplitudes develop. 
Most aerospace structures operate in random environment during their service life, i.e. 
the applied excitation is not deterministic. The following are some of the examples for 
random excitation. 

• The gust excitation caused by atmospheric turbulence 

• The excitation due to jet and rocket noise 

• The aerodynamic excitation due to boundary layer turbulence and fluctuating 
wake forces. 

• The ground loads induced during taxi, takeoff and landing of an aircraft or 
during the transportation of a launch vehicle to the launch pad. 

These results in need of modeling the system properties and excitation as random. 

In the present study, the system properties like longitudinal modulus, transverse modulus 
and Poison’s ratio, as well as the excitation are considered to be random. 
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In this section, a detailed study is presented to obtain the second order statistics of the 
response for sandwich honeycombed cylindrical, spherical, and conical panels using 
stochastic classical approach adopting the first order perturbation theory. 

3.4.2 FORMULATION: 


The displacements satisfying all the boundary conditions can be expressed as 
uix,y,t) = cos^^sin^exp(iu)0 


/n=l n=l 

CO 00 


vix, y, 0 = E E — cos exp(z cur) 

a b 


m-\ 

00 00 


w(x,y,t) = sin-^^sin^^exp(/(ur) 

a b 


m=l 


The transverse dynamic load in general can also be written as 


, . . mnx . .. . 

qix,y,t) = E sm sm-f-exp(z cur) 

a 0 

Substituting the above displacement forms in (2.25), results in 


^11 ^12 ^13 

^21 ^22 ^23 

L^31 ^32 ^33, 


■fA,' 


qx 

A, 

► zr * 

qi 

.k. 


.^3. 


(338) 


(3.39) 


All the random quantities are expressed as 

a,j = + aj and q. = + q! where i, j=l,2 and 3. 

Substituting all the random quantities in the equation (3.39), and collecting the zero order 
terms, we get 


d d d 

^11 ^12 ^13 

d d d 

^21 ^22 ^23 

d 


^31 ^32 


*33 


/ — — 


r 


» =z < 

d 

qi f 

.k. 


d 

93 

V J 


(3.40) 


Solving the above deterministic equation (3.40) we can get the deterministic response. 
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Collecting the first order terms, we get 


d d d\\.r 

^11 ^12 ^13 r^i 


d d 

I ^23 


^33 ^3'' 


A r r r r 1 d 

1 ^11 ^12 ^13 

A r , r r r a d 

2 ?“+ ^21 ^22 ^23 j ‘^2 

A r r r r f, d 

3 *^31 ^32 ^3 


(3.41) 


Using Taylor’s series expansion 


da. 

• 'J u '■ 




r SA," , 

a;=z-^; + 


/ 36/ * a?/ 


(3.42) 


/=1,2 and 3, i, j=l,2 and 3 


Substituting the equation (3.42) in (3.41), the unknowns are — ^ and — By 


solving equations (2.36) and (2.38), we can get the above said unknowns. Therefore, die 
total response becomes 


Aj=Aj + p, + T-T 

‘ obi * 


(3.43) 


The variance of the response can be evaluated as follows 


Var(A,)= 4 A;-A/| 


Var(A;)=£ + 

I / dbi dq, 


= P^ctEi;^ + +5V?/' 


(3.44) 


dA. 

where P = — - 


“"Y f^A/Y f^A/Y f^A/Y 

T and5=-;^ 


5D(A0 = +!2'o£22'' 


(3.45) 


i.e. Standard deviation of the response is the square root of variance of the response. 
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3.5 ANALYSIS OF HONEYCOMB 


This section deals with the evaluation of equivalent material properties of honeycomb for 
stability and dynamic analysis of laminated sandwich panels. The In-plane equivalent 
core material properties of honeycomb have been defined as [23] 

The equivalent modulus parallel to x, 

\2EJa^ cosO 


E, =■ 


a^{h-\-l sin 9)h!' sin * 6 


The equivalent modulus parallel to y, 


^ _ 1 2EJan (/? + /sin i9) 
' cr, W’ cos^ 0 


(3.46) 


(3.47) 


The Poisson’s ratio’s are calculated by 


-> — 


COS' (9 


h 

(-j + sind)sm& 


= 


h 

(y +sin^)sin^ 


COS' 9 


(3.48) 


The shear modulus can be obtained as G,, = 


UEJih+lsina) 


Ibh' cos9(l + 2h) 

For a cell with uniform thickness ‘t’, the above properties reduces to 




E, oos6 


^ (y + sin^)sin'^ 


h 


E^{-+sm9) 
_ {L\i I 


E,-(-) 


I cos 9 


(3.49) 


E,a + sin9) 


r 2h. 

(y)-COS^(l + y) 


(3.50) 


where Ei= Longitudinal modulus of honeycomb core, 
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E 2 = Transverse modulus of honeycomb core, 

Es= Modulus of elasticity of the solid from which the cell walls are made, 
Vy = Poisson’s ratios, 

/ = Inclined length of cell wall, 
t = Cell wall thickness, 

0 = Cell wall angle, 
h = Height of side wall of cell 
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CHAPTER 4 


RESULTS AND DISCUSSIONS 

4.1 PRESENT WORK 

Results have been obtained for cylindrical, spherical and conical shapes as particular cases 
of the general formulation in the previous chapter. Honeycomb sandwich laminated shell 
dimensions and lay up sequence are as shown in figure 4.1. Composite laminate is 
assumed to be axisymmetric with respect to shell mid surface. “B” matrix is zero due to 
symmetry. Moreover, the shell is assumed to be undergoing axisymmetric deformations, 
so the displacements and rotations are independent of the circumferential coordinate. 
Hence Vo drops from the formulation. 

Mathematically this assumption of axisymmetric displacement can again be written 

a 



Figure 4.1 : Lay up of the composite laminate (all the dimeiaions are in ‘mm’) 
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Substituting the above conditions (4.1) in the equations of motions of cylinder, 
sphere and conical shells the system coefficients take the following form for the specific 
shapes. 

4.1.1 CYLINDRICAL SHELLS: 


m 7t . 
— 

a 


a, 


niTt 

aR 


^12 


mn 


a 

33 ^2 


(4.2) 


4.1.2 SPHERICAL SHELLS: 


*^11 ~ 2 ^11 
a 


~ n (^11 -^12 ) 

aR 


, A A \ 

aR 


(/I,, + 2/4,, + /4„ ) 

«33 = ^ 


R- 


( 4 . 3 ) 


4.1.3 CONICAL SHELLS: 


a,, — ■ 


2 2 
m n 


4, 


^^13=- 


mjv 

aR, 


^12 


^31 =- 


rriK 

aR, 


^12 


■^22 

^33 =-^ 

r} 


(4.4) 


Where Rc is the radius of the cone. It is a function of ‘x’ i.e. Rc = R (x). 

The sandwich shell is modeled as an equivalent laminated construction as discussed 
in chapter 3. The assumed mean values of input variables used for graphite-epoxy and 
equivalent properties of aluminum core are given in Table 1. Table 2 represents the 
dimensions of the shells. 
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Table 1: Characteristics of primary variables 


Table 2: Dimensions of shells 



Length (m) 

Radius (m) 

Cylinder 

3.624 

2.000 


1.176 

jjHffljll 

Cone 1 

2.920 

2.00-1.176 

Cone 2 

0.980 

2.00-1.400 


Input 

Variable 

Mean values of 
Graphite-epoxy 

Mean values of 
Aluminum 

En 

3.05810 Gpa 

0.00466 Gpa 

E 22 

0.06191 Gpa 

111111 ^ 111 ^ 11 ^ 11 ^^ 

^2 

0.346 

0.330 

Gi2 

0.05088 Gpa 

0.00129 Gpa 

P 

1.7gm/cc 

0.0336gin/cc 


4.2.BUCKLING DUE TO AXIAL COMPRESSION OF SHELLS: 

The mean and the variance of the critical buckling load of cylinder, sphere and 
conical shells can be detemined by solving equations (3.5) and (3.19), using equations 
(4.2), (4.3) and (4.4) respectively. Here ‘NxCr’ is the Critical buckling load in ‘x’ direction 
and ‘a’ is the length of the shell. 

Numerical results are evaluated for the critical compressive load per unit 
circumference Nx for different types of shells considered 
4.2.1 THE MEAN BUCKLING LOADS 

Mean critical buckling load is determined from the equation (3.5) using equations 
(4.2), (4.3) and (4.4) for cylinder, sphere and conical shells respectively. The results for the 
mean of the critical buckling load (Nx) are presented in table 3 
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Table 3: Mean critical buckling loads for different shells 


Shell 

Critical buckling load (N/m) 

Cylinder 

318908.78125 

Sphere 

97129.26562 

Conel 

207040.75000 

Cone2 

23320.73829 


4.2.2 SD OF BUCKLING LOADS 

Equation (3.20) gives the Standard Deviation of critical buckling load, Nx as a 
function of standard deviation of input random variables. The ratios of SD to mean for the 
material properties are taken as 0, 5, 10, 15 and 20 percent. The three inputs are varied 
independently and the results are presented five sets in tabular form. 

Table 4 represents the change of SD/Mean of cylinder buckling load in five sub 
tables. Each table has one variation coefficient for E 22 variation 

coefficients assuming five values along the rows and the columns, respectively. Table 5 
and 6 represent the corresponding values for the spherical and the conical shells. 

It is inferred firom the results that 

• For the cone, the buckling load is position dependent as the radius changes along its 
length. The buckling load is observed to be increasing with decreasing radius of 
the shell due to increase of stiffixess. 

• All the shells are more strongly affected with changes in transverse modulus (E 22 ) 
than longitudinal modulus (En). 
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Table 4: SD/Mean of cylinder buckling load to SD/Mean of all RV’s 

sd/mean of y ,2 =0.0 sd/mean of u,, =0.05 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00007 

0.01536 

0.03073 

0.04609 

0.06146 

0.05000 

0.06146 

0.06769 

0.07273 

0.08044 

0.09013 













0.20000 

0.2637 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 







0.05000 

0.06592 

0.06769 

0.07273 

0.08044 


0.10000 

0.13185 

0.13274 

0.13538 

0.13967 

0.14547 







0.20000 

0.2637 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean of y,, =0.10 sd/mean of y,, =0.15 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.2CK>00 

0.00000 

0.00015 

0.01536 

0.03073 

0.04609 

0.06146 

0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 

0.13185 

0.13274 

0.13538 

0.13967^ 

0.14547 






0.2071 








sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00022 

0.01537 

0.03073 

0.04609 

0.06146 

0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 

0.13185 

0.13274 

0.13538 

0.13967 

0.14547 

0.15000 

0.19777 

0.19837 




0.20000 

0.2637 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00030 





0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 

0.13185 

0.13274 

0.13538 

0.13967 

0.14547 

0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.2071 

0.20000 

0.26370 

0.26415 

0.26548 

0.2677 

0.27077 


sd/mean of y ,2 = 0.20 
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Table 5: SD/Mean of sphere buckling load to SD/Mean of all RV’s 


sd/mean of t >,2 = 0.0 


sd/mean of t;,, = 0.05 


sd/mean 


■■ 















HQ 

0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 

0.13185 

0.13185 

0.13538 

0.13967 


0.15000 

0.19777 

0.19837 

0.20015 

0.20307 


0.20000 

0.26370 

0.26415 

0.26548 

0.26770 

0.27076 




■■ 










0.00000 

0.00007 

0.01536 

0.03073 

0.04609 

0.06146 



0.06769 

0.07273 

0.08044 

0.0901 S 


0.13185 

0.13185 

0.13538 

0.13967 

0.14547 






0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.26770 

0.27076 


sd/mean of i),^ =0.10 


sd/mean of y,, =0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00015 

0.01537 

0.03073 

0.04609 

0.06146 

0.05000 

0.06592 

I 

0.06769 

0.07273 

0.08044 







QQ9 

0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.20710 








sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00022 

0.01537 

0.03073 

0.04609 

0.06146 

0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 


■ 



0.14547 

0.15000' 





0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.26770 

0.2707^ 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00030 

0.01537 

0.03073 

0.04609 

0.06146 

0.05000 

0.06593 

0.06769 

0.07274 

0.08044 

0.09013 





0.13967 

0.14547 

0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.26707 

0.27077 


sd/mean of = 0.20 
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Table 6: SD/Mean of cone buckling load to SD/Mean of all RV’s 


sd/mean of = 0.0 


sd/mean of Uj, = 0.05 













0.06146 

0.05000 

0.06592 

0.06769 

1 

0.07273 

0.08044 


0.10000 

0.13185 

0.13185 

0.13538 

0.13967 


0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.2000( 

0.00000 







0.06592 

0.06769 

0.07273 

0.08044 

0.09015 





0.13967 

0.14547 





0.20307 

0.2071C 


0.26370 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean of l >,2 =0.10 


sd/mean of u,, =0.15 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.200(X) 

0.00000 

0.00015 

0.01537 

0.03073 

0.04609 

0.06146 

0.05000 

0.06592 

0.06769 

0.07273 

0.08044 

0.09013 

0.10000 

0.13185 

0.13185 

0.13538' 

0.13967 

0.14547 

0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.2677 

0.27076 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00022 

0.01537 

0.03073 

0.04609 

0.06146 







0.10000 

0.13185 

0.13274 

0.13538 

0.13967 

0.14547 






0.2071 Q 

0.20000 

0.26370 

0.26415 

0.26548 

0.2677 

0.2707?! 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.0003 

0.01537 

0.03073 

0.04609 

0.06146 













0.15000 

0.19777 

0.19837 

0.20015 

0.20307 

0.20710 

0.20000 

0.26370 

0.26415 

0.26548 

0.2677 

0.27077 


sd/mean of y ,2 = 0.20 
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• All the shells are least affected with changes in the Poisons ratio 0^2 - 

• The SD/Mean of buckling load is linearly varying with the SD/Mean of the input 
random variables considered. 

4.3. FREE VIBRATION OF SHELLS: 

The natural frequencies of cylinder, sphere and conical shells have been determined by 
setting the external loads and in-plane loads as zeros in equation (2.10) and by solving 
equation (3.27) using (4.2), (4.3) and (4.4) respectively. The Variance and Standard 
deviation of o^have been determined by equations (3.36) and (3.37). As mentioned 
earlier the deflection is assumed to be axisymmetric. 

4.3.1 MEAN NATURAL FREQUENCIES 

In this case, the characteristic equation, (equation (3.27)), for the mean of natural 
frequency is quadratic. Hence, there are two roots corresponding to the value of X . 
Table7 gives the mean of first two natural frequencies of cylinder, sphere and conical 
shells corresponding to first five wave numbers.. 


Table?: Mean of Natural Frequencies (HZ) 


Mode 

Frequency 

Cylinder 

Sphere 

Conel 

Cone2 

1 


335.7443 

535.4978 

339.02500 

M3.13251 


2^ 



790.5953 

2327.45337 

2 

1“ 

342.0127 

571.5809 

342.5817 

343.45821 


2^9 

1262.898 

3960.036 

1564.7740 

4650.49658 

3 

1“ 

342.9012 

578.6065 

343.1382 

343.51700 


2*^ 

1889.439 



6974.54346 

4 

i« 

343.1963 



343.53992 


2^ki 

2517.086 



9298.83203 

5 

1“ 

343.3300 

582.2404 

343.4131 

343.54565 


2 ^ 

3145.131 

9718.854 

3902.469 

11623.21680 
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4.3.2 SD OF NATURAL FREQUENCIES 


The SD of the natural frequencies is obtained as a function of SD of input random 
variables from the equation (3.37). The ratio of SD to mean for the material properties is 
varied as mentioned earlier for the buckling study. 

The Table (8) to table (22) represents the change of SD/Mean of square of natural 
frequencies with respect to the SD/Mean of random variables for all shells up to mode 5. 
The normalized results for conel and cone2 are identical. The following conclusions can 
be made from the results 

• Cylinder and spherical shells are more strongly affected with changes in 
longitudinal modulus (Ei i) than transverse modulus (E 22 )- 

• Conical shells show greater sensitivity with changes in transverse modulus (E 22 ) 
than longitudinal modulus (En). 

• All the shells are least affected with changes in the Poisons ratio t>, 2 . 

• As the mode number increases, the influence of Eu is decreasing and that of E 22 is 
increasing for all the shells. 

• The SD/Mean of square of the natural frequency for conel and cone2 is same due 
to the normalization scheme adopted. 

• The SD/Mean of square of the natural frequency increases linearly with the input 
SD/Mean of random variables for all frie shells. 

• The stiffiiess of conical shell is position dependent as the radius changes along its 
length. The natural frequencies are observed to be increasing with decreasing 
radius of the shell. 

• The natural frequencies increase with increase in the mode number. 



Table (8): SD/Mean of square of natural frequency of cylinder for model to 

SD/MeanofallRV’s 


sd/mean ofu ,2 = 0.05 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00029 

0.20598 

0.41195 

0.61793 

0.82391 

0.05000 

0.06761 

0.21679 

0.41746 

0.62162 

0.82668 

0.10000 

0.13521 

0.24639 

0.43557 

0.63255 

0.83493 

0.15000 

0.20282 

0.28907 

0.45917 

0.65036 

0.84850 

0.20000 

0.27042 

0.33993 

0.49278 

0.67451 

0.86715 


sd/mean =0.10 


Sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00058 



0.20598 

0.41195 

0.61793 

0.82391 






0.82668 

0.10000 

0.13521 

0.24639 

0.43557 

0.63255 

0.83493 

0.15000 

0.20282 

0.28907 

0.45917 

0.65036 

0.84850 

0.20000 

0.27042 

0.33993 

0.49278 

0.67451 

0.86715 


sd/mean of t >,2 = 0. 15 











0.10000 

0.15000 

0.20000 






0.82391 

0.05000 

0.06761 

0.21679 

0.41746 

0.62162 

0.82668 





0.63255 

0.83493 











0.67451 

0.86715 
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Table (9): SD/Mean of square of natural frequency of cylinder for model to 


SD/Mean of all RV’s 


sd/mean ofy ,2 = 0.05 


Sd/mean 


i 




0.00000 

0.05000|0. 10000 

0.15000 

0.20000 

0.00000 

0.00014 

0.01170|0.02340 

0.03510 

0.04680 

0.05000 

0.07276 

0.073700.07643 

0.08079 

0.08651 

0.10000 

0.14552 

0.14599! o . 14739 

! 

0.14970 

0.15286 

0.15000 

0.21829 

0.218600.21954 

1 

0.22109 

0.22325 

0.20000 

0.29105 

0.291280.29199 

0.29316 

0.29479 


sd/mean ofy,, =0.10 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00029 

0.01170 

0.02340 

0.03510 

0.04680 

0.05000 

0.07276 

0.07370 

0.07643 

0.08079 

0.08651 

0.10000 

0.14552 

0.14599 

0.14739 

0.14970 

0.15286 

0.15000 

0.21829 

0.21860 

0.21954 

0.22109 

0.22325 

O 

o 

o 

o 

C\| 

o 

0.29105 

0.29128 

0.29199 

0.29316 

0.29479 


sd/mean ofy,, =0.15 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 













0.10000 

0.14552 

0.14599 

0.14739 

0.14970 

0.15286 

0.15000 

0.21829 

0.21860 

0.21954 

0.22109 

0.22325 

0.20000 

0.29105 

0.29128 

0.29199 

0.29316 

0.29479 
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Table (10): SD/Mean of square of natural frequency of cylinder for mode 3 

to SD/Mean of all RV’s 

sd/mean ofu ,2 = 0.05 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00021 

0.03996 

0.07992 

0.11988 

0.15984 

0.05000 

0.07344 

0.08361 

0.10854 

0.14059 

0.17590 

0.10000 

0.14687 

0.15221 

0.16721 

0.18959 

0.21707 

0.15000 

0.22031 

0.22391 

0.23436 

0.25082 

0.27219 

0.20000 

0.29375 

0.29645 

0.30443 

0.31727 

0.33442 


sd/mean of = 0.10 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00042 

0.03996 

0.07992 

0.11988 

0.15984 

0.05000 

0.07344 

0.08361 

0.10854 

0.14059 

0.17590 







0.15000 

0.22031 

0.22391 

0.23436 

0.25082 

0.27219 

0.20000 

0.29375 

0.29645 

0.30443 

0.31727 

0.33442 


sd/mean ofv ^2 = 0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00063 

0.03996 

0.07992 

0.11988 

0.15984 










0.16721 








0.27219 

0.20000 

0.29375 

0.29645 

0.30443 

0.31727 

0.33442 
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Table (11): SD/Mean of square of natural frequency of cylinder for mode 4 

toSD/MeanofallRV’s 


sd/mean oft >,2 = 0.05 


E22 ^ 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00023 

0.04918 

0.09836 

0.14753 

0.19671 

0.05000 

0.07366 

0.08857 

0.12288 

0.16490 

0.21005 

0.10000 

0.14732 

0.15531 

0.17713 

0.20849 

0.24756 

p 

01 

0 

0 

0 

0.22097 

0.22638 

0.24187 

0.26570 

0.29585 

0.20000 

0.29463 

0.29871 

0.31062 

0.32951 

0.35427 


sd/mean of u ,2 =0.10 







■■ 







0.00000 

0.00047 

0.04918 

0.09836 

0.14753 

0.19671 

0.05000 

0.07366 

0.08857 

0.12288 

0.16490 

0.21005 

0.10000 

0.14732 

0.15531 




0.15000 

0.22097 

0.22638 



0.29585 

0.20000 

0.29463 

0.29871 





sd/mean oft >,2 = 0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 






0.19671 

0.05000 

0.07366 

0.08857 

0.12288 

0.16490 

0.21005 

0.10000 

0.14732 

0.15531 

0.17713 

0.20849 

0.24756 

0.15000 

0.22097 

0.22638 



0.29585 

0.20000 

0.29463 

0.29871 



0.35427 


\ 
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Table (12): SD/Mean of square of natural frequency of cylinder for mode 5 

to SD/Mean of all RV’s 


sd/mean oft >,2 = 0.05 


"si^ean 
E22 ^ 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 






0.05000 

0.07376 

0.09102 

0.12967 

0.17617 

0.22570 





0.21761 

0.25935 

nn^ 


HQ 




0.20000 

0.29503 

0.29891 

0.31372 

0.33562 

0.36407 


sd/mean of y,, = 0. 10 


sd/mean 







0.00000 

0.050)0 

0.10000 

0.15000 

0.20000 

0.00000 

0.00049 

0.05333 

0.10666 



0.05000 

0.07376 

0.09102 

0.12967 

0.17617 

0.22570 




0.18203 




0001 





m 

0^00 






sd/mean ofyj 2 =0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00073 

0.05333 

0.10666 

0.15998 

0.21331 







0.10000 

0.14752 

0.15686 

0.18203 

0.21761 

0.25935 

p 

cn 

0 

0 

0 

0.22127 

0.22761 



0.30735 

0.20000 

0.29503 

0.29891 

0.31372 

0.33562 

0.36407 
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Table (13): SD/Mean of square of natural frequency of sphere for mode 1 

to SD/Mean of all RV’s 


sd/mean oft >,2 = 0.05 



0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

■g 





1.19572 





0.89879 

1.19722 

0.10000 

0.11966 

0.32199 

0.60972 

0.90474 

1.20170 

0.15000 

0.17948 

0.34867 

0.62422 

0.91458 

1.20912 

0.20000 

0.23931 

0.38292 

0.64398 

0.92817 

1.21943 


sd/mean of D,, = 0. 10 


Sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00373 

0.29895 

0.58789 

0.89880 

1.19573 

0.05000 

0.05994 


0.60086 

0.89879 

1.19722 

0.10000 

0.11971 

0.32201 

0.60973 

0.90475 

1.20170 

0.15000 

0.17951 

0.34869 

0.62423 

0.91458 

1.20912 

0.20000 

0.23933 

0.38293 

0.64398 

0.92818 

1.21944 


sd/mean ofv^^ - 0-15 



|H|| 








0.10000 

0.15000 

0.20000 






1.19574 

0.05000 

0.06009 

0.30491 

0.60087 

0.89880 

1.19723 

0.10000 

0.11978 

0.32204 

0.60974 

0.90476 

1.20171 

0.15000 

0.17956 

0.34871 

0.62424 

0.91459 

1.20913 

0.20000 

0.23936 

0.38296 

0.64400 

0.92819 

1.21945 
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Table (14); SD/Mean of square of natural frequency of sphere for mode 2 

toSD/MeanofallRV’s 


sd/mean of v ^2 = 0.05 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00178 

0.08515 

0.17026 

0.25339 

0.34051 

0.05000 

0.06440 

0.10674 

0.18203 

0.26337 

0.34654 

0.10000 

0.12876 

0.15435 

0.21346 

0.28600 

0.36404 

0.15000 

0.19313 

0.21106 

0.25746 

0.32018 

0.39146 

0.20000 

0.25750 

0.27120 

0.30869 

0.36266 

0.42691 


sd/mean ofL >,2 =0.10 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00356 

0.08520 

0.17029 



0.25541 

0.34053 











0.28602 

0.36405 

0.15000 

.019315 

0.21108 

0.25748 

0.32020 

0.39148 


0 ^^ 


0.30871 

0.36268 

0.42692 


sd/mean oft >,2 =0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00535 

0.08529 

0.17034 






0.18210 

0.26342 

0.34658 




0.21352 

0.28605 

0.36407 




0.25751 



B 



0.30873 

0.36270 

0.42694 
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Table (15): SD/Mean of square of natural frequency of sphere for mode 3 

to SD/Mean of all RV’s 


sd/mean ofL >,2 = 0.05 


sd/mean 

0.00000 

0.05000 

O.ICKKK) 

0.15000 

0.20000 

0.00000 

0.00176 

0.04614 

0.09222 

0.13832 

0.18442 

0.05000 

0.06523 

0.07998 

0.11295 

0.15292 

0.19561 

0.10000 

0.13043 

0.13834 

0.15973 

0.19011 

0.22588 

0.15000 

0.19564 

0.20100 

0.21628 

0.23959 

0.26886 

0.20000 

0.26085 

0.26489 

0.27667 

0.29525 

0.31945 


sd/mean of = 0. 10 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 







0.05000 

0.06531 

0.07994 

0.11299 

0.15295 

0.19564 

0.10000 

0.13047 

0.13837 

0.15976 

0.19014 

0.22590 







0.20000 

0.26087 

0.26491 

0.27668 

0.29526 

0.31947 


sd/mean ofi >,2 =0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00529 

0.04641 

0.09236 

0.13841 

0.18449 

0.05000 

0.06542 

0.08004 

0.11306 

0.15300 

0.19568 





0.19018? 

0.22594 

0.15000 

0.19570 

0.20106 

0.21634 

0.23964 

0.26890 

0.20000 

0.26089 

0.26494 

0.27671 

0.29529 

0.31949 
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Table (16): SD/Mean of square of natural frequency of sphere for mode 4 

to SD/Mean of all RV’s 


sd/mean of = *^-05 


sd/mean 






^22\^E^ 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 






Qlllllllll 

0.05000 

0.06550 

0.07318 

0.09241 

0.11767 

0.14587 

0.10000 

0.13101 

0.13500 

0.14632 

0.16346 

0.18479 

0.15000 

0.19653 

0.19919 

0.20703 

0.21948 

0.23580 

0.20000 

0.26201 

0.26403 

0.26999 

0.27965 

0.29263 


sd/mean of u,, = 0. 10 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

. 

0.20000 

0.00000 

0.00351 

0.03277 

0.06526 

0.09781 

0.13037 

0.05000 

0.06552 

0.07324 

0.09246 

0.11771 

0.14590 

0.10000 

0.13105 

0.13504 

0.14635 

0.16348 

0.18482 







0.20000 

0.26203 

0.26404 

0.27001 

0.27966 

0.29265 


sd/mean ofy,, =0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 




0.06537 

0.09788 

0.13043 

0.05000 

0.06559 

0.07335 




0.10000 

0.13111 

0.13509 


0.16353 

0.18486 




0.20709 

0.21953 

0.23585 

0.20000 

0.26205 

0.26407 

0.27003 

0.27969 

0.29267 
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Table (17): SD/Mean of square of natural frequency of sphere for mode 5 

to SD/Mean of all RV’s 


sd/mean oft ), 2 = 0.05 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00175 

0.02643 

0.05278 

0.07915 

0.10552 







0.10000 

0.13128 

0.13390 

0.14148 

0.15328 

0.16842 

0.15000 

0.19691 

0.19867 

0.20386 

0.21222 

0.22340 

0.20000 

0.26254 

0.26387 

0.26779 

0.27421 

0.28295 


sd/mean of d ,2 = 0. 10 


Sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 







0.05000 

0.06573 

0.07324 

0.08428 

0.10287 

0.12430 

0.10000 

0.13132 

0.13394 

0.14152 

0.15331 

0.16845 

0.15000 

0.19694 

0.19869 

0.20388 



0.20000 

0.26256 

0.26388 

0.26781 




sd/mean ofUjj = 0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 






0.10564 

0.05000 

0.06584 

0.07335 

0.08437 

0.10294 

0.12437 





0.15336 

0.16850 

0.15000 

0.19697 

0.19873 

0.20392 

0.21227 

0.22345 

0.20000 

0.26259 

0.26391 

0.26784 

0.27426 

0.28299 
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Table (19): SD/Mean of square of natural frequency of cone for mode 2 

to SD/Mean of ail RV’s 

sd/mean of?;, 2 = 0.05 


^^mean 






E22 tiY- 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00001 

0.05305 

0.10611 

0.15916 

0.21221 

0.05000 

0.06506 

0.08395 

0.12446 

0.17194 

0.22196 

0.10000 

0.13012 

0.14052 

0.16790 

0.20558 

0.24893 

0.15000 

0.19518 

0.20226 

0.22216 

0.25185 

0.28832 

0.20000 

0.26024 

0.26560 

0.28104 

0.30505 

0.33580 


sd/mean ofv,, =0.10 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00001 

0.05305 

0.10611 

0.15916 

0.21221 

0.05000 

0.06506 

0.08395 

0.12446 

0.17194 

0.22196 

0.10000 

0.13012 

0.14052 

0.16790 

0.20558 

0.24893 







0.20000 

0.26024 

0.26560 

0.28104 

0.30505 

0.33580 


sd/mean oft),, = 0.15 






■■ 








0.00000 

0.00002 

0.05305 

0.10611 



0.05000 

0.06506 

0.08395 

0.12446 

0.17194 

0.22196 












0.28832 

0.20000 

0.26024 

0.26560 

0.28104 

0.30505 

0.33580 
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Table (20): SD/Mean of square of natural frequency of cone for mode 3 to 

SD/Meanof allRV’s 


sd/raean ofy,, = 0.05 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00004 

0.03151 

0.06302 

0.09452 

0.12603 

0.05000 

0.06552 

0.07271 

0.09091 

0.11501 

0.14205 

0.10000 

0.13105 

0.13478 

0.14541 

0.16158 

0.18182 

0.15000 

0.19657 

0.19908 

0.20643 

0.21812 

0.23350 

0.20000 

0.26210 

0.26398 

0.26956 

0.27862 

0.29082 


sd/mean of w,, = 0. 1 0 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00008 

0.03151 

0.06302 

0.09452 

0.12603 

0.05000 

0.06552 

0.07271 

0.09091 

0.11501 

0.14205 


0.13105 


0.14541 



0.15000 

0.19657 

0.19908 

0.20643 

0.21812 

0.23350 

O 

§ 

d 

0.26210 

0.26398 

0.26956 

0.27862 

0.29082 


sd/mean ofv ,2 = 0.15 


sd/mean 







0.00000 

0.05000 

0.10000 

0.15000 

0.20000 












0.14205 


0.13105 

0.13478 

0.14541 

0.16158 

0.18182 

0.15000 

0.19657 

0.19908 

0.20643 

0.21812 

0.23350 

0.20000 

0.26210 
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Table (21); SD/Mean of square of natural frequency of cone for mode 4 to 

SD/MeanofallRV’s 


sd/meah ofu ,2 = 0.05 


Sd/mean 

' e 22 \.^ E ^ 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00005 

0.02433 

0.04865 

0.07298 

0.09730 

0.05000 

0.06568 

0.07004 

0.08173 

0.09818 

0.11739 

0.10000 

0.13136 

0.13359 

0.14008 

0.15027 

0.16347 







0.20000 

0.26271 

0.26384 

0.26718 

0.27266 

0.28015 


sd/mean ofu,, =0.10 


Sd/mean 







0.00000 

0.05000 

0-10000 

0.15000 

0.20000 

0.00000 

0.00011 

0.02433 

0.04865 

0.07298 

0.09730 







0.10000 

0.13136 

0.13359 

0.14008 

0.15027 

0.16347 

0.15000 

0.19704 

0.19853 

0.20295 

0.21012 

0.21975 

0.20000 

0.26271 

0.26384 

0.26718 

0.27266 

0.28015 


sd/mean = 0.15 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 

0.00000 

0.00016 

0.02433 

0.04865 

0.07298 

0.09730 






0.11739 







0.15000 

0.19704 

0.19853 

0.20295 

0.21012 

0.21975 


0 ^^ 
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Table (22); SD/Mean of square of natural frequency of cone for mode 5 to 

SD/Mean of all RV’s 


sd/mean ofi >,2 = 0.05 


■ 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 




0.04213 

0.06319 

0.08426 

0.05000 

0.06575 

0.06904 

0.07809 

0.09119 

0.10688 

0.10000 

0.13150 

0.13317 

0.13808 

0.14589 

0.15618 







0.20000 

0.26299 

0.26384 

0.26635 

0.27048 

0.27616 


sd/mean of D,, = 0.10 





■■ 









0.00000 

0.00012 

0.02106 

0.04213 

0.06319 

0.08426 

0.05000 

0.06575 

0.06904 

0.07809 

0.09119 

0.10688 

0.10000 

0.13150 

0.13317 

0.13808 

0.14589 

0.15618 






0.21449 

0.20000 

0.26299 



0.27048 

0.27616 


sd/mean ofw ,2 = 0.15 


sd/mean 

0.00000 

0.05000 

0.10000 

0.15000 

0.20000 






0.08426 

0.05000 

0.06575 

0.06904 

0.07809 

0.09119 

0.10688 

0.10000 






0.15000 

0.19725 

0.19837 

0.20170 

0.20712 

0.21449 

0.20000 

0.26299 

0.26384 

0.26635 

0.27048 

0.27616 
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4.4 FORCED VIBR4TI0N OF SHELLS: 


The deflection amplitudes of cylinder, sphere and conical shells have been 
determined by solving the equation (3.40). The Variance and Standard deviation of 
deflections are determined from the equations (3.44) and (3.45). 

Numerical results are presented for the forced vibrations of cylinder, sphere and 
conical shells. The boundary conditions for the shell are considered to be simply 
supported. The deflections, as mentioned earlier are assumed to be axisymmetric. All the 
SD of response is normalized with their respective mean values. The assumed mean values 
of the input primary variables are given in Tablel. The mean external load is taken as 
5000 N/sq.mm. 

4.4.1 MEAN RESPONSE AMPLITUDE 

The mean response can be determined by solving equation (3.40). Table (23) gives 
the mean response of cylinder, sphere and conical shells. 

Table (23): Mean Deflection amplitude 


Shell 

U - Deflection (m) 

W - Deflection (m) 

Cylinder 

0.00387 

0.02103 

Sphere 

0.00303 

0.00723 

Conet 

0.00311 

0.02103 

Cone2 

0.00104 

0.02102 


4.4.2 SD OF RESPONSE 

The SD of the deflection amplitude is obtained as a function of SD of input random 
variables from the equation (3.45). The ratio of SD to mean for the material properties and 
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the excitation are considered to chano<» j , 

to change mdepcndently. The results are obtained for their 

values as 0, 5, 10, 15 and 20 percent. 

The Figures 4.2 to 4.5 represent the change of SD/Mean of response with respect to 
the SD/Mean of input random variables. 


It is observed from the results that 

. AU the sheUs are more strongly affected with chtmges in longitudinal modulus (E„) 

and the external loading (Q) compared to changes in transverse modulus (E.) for 
deflection - U 

All the shelU are more sensitive to changes in transverse modulus (E,,) and 
external loading (Q) than the longitudinal modulus (E„) for deflection - W. 

The mean response is decreasing with decreasing of radius of conical shell due to 
change in stiffiiess along its length. 

• All the shells are least affected with changes in the Poison’s ratio 0 , 2 . 

• The SD/Mean of deflection for conel and cone2 is same due to normalization. 

The SD/Mean of response increases linearly with the input SD/Mean of random 
variables for all the shells. 

• The SD/Mean of deflections increases with increase in mode number for all the 
shells. 
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sd/mean of deflection (u) sd/mean of deflection (u) 


3.5 



0 0.05 0.1 0.15 0.2 0.25 

sd/mean of all RV’s 



sd/mean of all RV’s 

Figure 4.3: Variation of sd/mean of deflections for spherical panel 
(a) u - deflection (b) w - deflection 
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sd/mean of all RV’s 

Figure 4.4; Variation of sd/mean of conical panel-1 
(a) u - deflection (b) w - deflection 
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sd/mean of deflection (w) sd/mean of deflection (u) 




sd/mean of all RV’s 

Figure 4.5: Variation of sd/mean of deflections for conical panel 2 
(a) u - deflection (b) w - deflection 
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4.5 SATELLITE PAYLOAD FAIRING: 


A typical satellite payload fairing considered here consists of 3 types of shells 
namely cylinder, sphere and cone of different dimensions (Figure 4.6). Results for 
independent stability and dynamic analysis of various shells have been presented in this 
chapter. But in the complete structure, these component shells are attached together and 
support each other. 

For the satellite payload fairing, two extremes of ideal support conditions may be 
considered for the shells to be simply supported and fixed ends. Simply supported 
condition allows rotation but not deflection but the fixed condition neither allows rotation 
nor deflection. The satellite payload fairing shell components will actually have flexible 
support conditions such that they allow some rotation and also some deflection. 

It can be conclude firom the results that, 

• The critical buckling load of the entire structure for the assumed support conditions 
is the lowest among all the critical buckling loads of shell gjven in table 3, which is 
for cone2. 

• The lowest natural fi'equency of the entire system for the assumed support 
conditions is the lowest among all the frequencies presented in the table 7, which is 
for cylinder. 

• The maximum deflections along x and z-axis respectively are the largest values 
among all the deflections presented in the table 23, which is for cylinder. 

• The critical buckling load, natural frequency and the deflections depends on the 
support conditions, dimensions of the shell, material properties and external load 
also. 


73 



8700 


For the actual conditions, the behavior of satellite payload fairing will be in 
between the behavior of the structure when the support conditions are simply 
supported and fixed. 
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Figure 4.6: Typical satellite payload fairing 


CHAPTER 5 


CONCLUSIONS AND SCOPE FOR FUTURE WORK 

5.1 CONCLUSIONS 

A general approach has been presented for stability and d5mamic analysis of 
laminated sandwich shells with random material properties and external loading. The 
second order statistics, that is, mean and standard deviation of buckling load, natural 
frequencies and forced response amplitude characteristics have been obtained for 
cylinder, sphere and conical shells as a frinction of known statistics of input random 
variables using perturbation technique. The method has been utilized to obtain numerical 
results for a specific shell configuration. In the buckling analysis, the buckling of shells 
under axial compressive load is studied. Axisymmetric vibration modes are considered in 
the vibrational analysis. Some conclusions that are based on this study are as follows; 

• SD/Mean of buckling load, natural frequency and deflection changes linearly with 
SD/Mean of input RV’s. 

• Larger response SD is associated with higher mean response. Hence for higher 
magnitude of loads, random analysis is more important for accurate evaluation of 
laminated shell response behaviour. 

• All the shells are least affected with changes in the Poisom ratio u,2 . 

• The buckling loads and natural frequencies show different sensitivity to 
randonmess in different basic material properties. The dominant material 
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property depends on the laminate construction, mode of deflection/vibration, and 
geometrical parameters. 

• The buckling load increases with decreasing radius of the shell. 

• All the shells in buckling are more strongly affected with changes in transverse 
modulus (E22) than longitudinal modulus (Eu). 

• For the free vibrations, Cylinder and spherical shells are strongly affected with 
changes in longitudinal modulus (Eu) and Conical shells show greater sensitivity 
with changes in transverse modulus (E22). 

• The effect of Ei 1 decreases and that of E22 increases with the increase in mode 
number in free vibration. 

• The ratio of SD/Mean of buckling load is same for cylinder, sphere and conical 
shell due to the normalization scheme adopted. 

• In forced vibration, the longitudinal deflection - u, all the shells are strongly 
affected with changes in En and Q and for radial deflection - w, all the shells are 
sensitive to the changes in E22 and Q. 

5.2 SCOPE FOR FUTURE WORK 

The present study is confined to analysis of thin shells (h/R< 0 . 1 ) so that classical 
thin shell theory is used. In this study, transverse shear deformation effects and rotatory 
inertia effects are neglected. The present study can be extended to analysis of thick shells 
including transverse shear deformation and rotatory inertia effects with different support 
conditions to develope a more complete understanding of composite shell behaviour. 
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The first order perturbation terms of random variables are only considered in the 
present study and higher order terms are neglected in comparison with their 
corresponding mean values. For future work, second and higher order terms can be 
considered to make the analysis more precise and valid in the nonlinear domain. 
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